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Environmental impact and depletion of mineral resources such as coal, oil and gas are prompting a reex¬ 
amination of an alternative to these resources. A safe and sustainable energy pathway which is crucial to 
sustainable development in addition to greenhouse gas emitters and its relationship with climate change 
are leading factors to look for adequate strategies concerning both energy saving and environmental pro¬ 
tection. Solar heat engines are attracting much interest nowadays and, as a consequence, different Stirling 
engine coupled to solar collector have been investigated since it meets the demands of the efficient use of 
energy and assuring environmental security. In recent years several prototypes and experimental facili¬ 
ties of solar Stirling engine have been developed. The future commercial interest of this alternative for 
electric power generation relies on a reduction of investment costs and on an increase of performance. 
The Stirling heat engine was first patented in 1816 by Robert Stirling. Since then, several Stirling engines 
based on his invention have been built in many forms and sizes. The engine works with a closed cycle and 
uses several gases as working fluid. Several prototypes have already been studied and produced but the 
alpha Stirling engine using the Ross Yoke linkage was not well studied although this kind of engine has 
the advantage of minimizing lateral forces acting on the pistons and leading to a more efficient and com¬ 
pact design compared to beta or gamma Stirling configuration, The objective of this work was the study of 
the effect of the geometrical and physical parameters on Ross Yoke Stirling engine performance in order 
to determine the significant thermodynamic parameters having an impact on engine performance. We 
proposed thermodynamic optimization of a Ross Yoke Stirling engine on the basis of a numerical model 
integrating the internal and external irreversibility. As a result, this analysis indicated that the perfor¬ 
mance of a Ross Yoke Stirling cycle engine with air as working gas depends critically on the geometrical 
parameters and heat input. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

Reduction of Greenhouse Gas (GHG) emissions can mainly be 
achieved through utilizing the abundant Renewable Energy 
Sources (RESs) and the implementation of Energy Efficiency (ENEF) 
measures. Renewable energy, with the availability of its renewabil- 
ity and non-pollution, will grow to be an effective and practical 
choice to guarantee the future development of the world. 

Developing renewable energy is an inevitable choice for sus¬ 
tainable economic growth, for the harmonious coexistence of hu¬ 
man and environment as well as for the sustainable 
development. Solar energy is one of the more attractive renewable 
energy sources, especially in MENA countries, and can be used as 
an input energy source for Stirling engine representing one of the 
most effective applications [1]. Stirling engine is a regenerative 
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externally heated engine operating with a cycle that has the same 
thermal efficiency with Carnot cycle. Since it can be powered by 
various heat sources and is also able to use solar radiation as en¬ 
ergy source. Theoretical and experimental numerous studies have 
been conducted on Stirling engines. Many types have been de¬ 
signed and built [2-10] however alpha Stirling engine using the 
Ross Yoke linkage was not well presented and studied although 
its advantage of minimizing lateral forces acting on the pistons 
and leading to a more efficient and compact design compared to 
others configuration. 

None of Stirling engine became competitive with the internal 
combustion engines dues to the irreversibility effects caused by 
external heating and imperfect regeneration. Several studies and 
prototypes were produced, in solar energy investigations. The first 
use of Stirling engine using solar energy was made by Parker and 
Malik in 1962. In their conversion system, solar radiation was con¬ 
centrated and reflected on to the heater of the engine by using a 
Freshnel lens. 
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Nomenclature 



A 

area (m 2 ) 

ch 

load 

Cp 

specific heat at constant pressure (J kg -1 K _1 ) 

cd 

conduction 

Cpr 

heat capacity of each cell matrix (W K _1 ) 

d 

expansion space 

Cv 

specific heat at constant volume (J kg -1 K _1 ) 

E 

entered 

d 

hydraulic diameter (m) 

ext 

outside 

D 

diameter (m) 

f 

cooler 

s 

regenerator efficiency 

h 

heater 

fr 

friction factor 

in 

irreversible 

J 

annular gap between displacer and cylinder (m) 

moy 

average 

G 

working gas mass flux (kg m -2 s _1 ) 

P 

loss 

k 

thermal conductivity (W m -1 K -1 ) 

Pa 

wall 

L 

length (m) 

pis 

piston 

M 

mass of working gas in the engine (kg) 

r 

regenerator 

m 

mass flow rate (kg s _1 ) 

rl 

regenerator cell 1 

m 

mass of gas in different component (kg) 

r2 

regenerator cell 2 

P 

pressure (Pa) 

S 

left 

Q 

heat (J) 

shtl 

shuttle 

Q 

power (W) 



R 

gas constant (J kg -1 1< -1 ) 

Greek letters 

T 

temperature (I<) 

6 

crank angle (rad) 

U 

convection heat transfer coefficient (W m -2 K _1 ) 

P 

working gas dynamic viscosity (kg m -1 s -1 ) 

V 

volume (m 3 ) 

P 

density (kg m -3 ) 

W 

work (J) 

CO 

angular frequency (rad s -1 ) 

Z 

displacer stroke (m) 

y 

Cp Cv 1 

Subscripts 



c 

compression space 




Popescu [11] studied the influence of the losses of heat and irre¬ 
versibility on the performance of the engine by finite-time thermo¬ 
dynamic optimization of an endo- exo-irreversible Stirling motor. 
He demonstrated that the most important reduction of the perfor¬ 
mances is due to the non-adiabatic regenerator. 

Kolin [12] studied several mechanisms of the piston movement 
transformation and showed that more the movement is discontin¬ 
uous more the real cycle approaches to the ideal cycle. 

In order to obtain a compact Stirling engine and with reduced 
price, Hirata [13] developed a prototype using the air as working 
fluid. He measured the mechanical losses for various values of load 
pressure to obtain the highest output. Karabulut and Koca [14] 
developed a Stirling engine with swept volume equal to 260 cm 3 . 
They studied the effect of load pressure and the heating tempera¬ 
ture on the torque. The maximum capacity of engine equal to 65 W 
was obtained for the heating temperature 1100°C and with the 
load pressure of 2, 5 bars. 

Formosa and Despesse [15] developed an analytical model for a 
free piston Stirling engine in order to investigate the effects of the 
technological and operating parameters on Stirling engine perfor¬ 
mance. He deduced that the cooler effectiveness affects the engine 
performances. 

Chin-Hsiang et al. [16] developed a numerical model for a beta- 
type Stirling engine with rhombic-drive mechanism by taking into 
account different thermal losses. They predicted all thermody¬ 
namic parameters. They concluded that the performance could 
be improved by adjusting the influential parameters including 
the regenerative gap, the distance between the two gears, the off¬ 
set distance from the crank to the center of the gear, and the heat 
source temperature. 

Karabulut [17] developed a numerical model to study a free pis¬ 
ton Stirling engine working with closed and open thermodynamic 
cycle. He concluded that the power exhibits very big variations 
with respect to the static position of the piston and the displacer, 


the hot end temperature, the piston damping coefficient, the dis¬ 
placer rod diameter and the stiffness of the piston spring. 

Cinar et al. [18] carried out an experimental study on a beta- 
type Stirling engine working at atmospheric pressure. He con¬ 
cluded that an increase in the hot-source temperature leads to 
an increase in the engine speed and torque and the power-output. 
Besides, the performance can be increased by using a working fluid 
with a higher thermal conductivity, such as helium or hydrogen. 

Cullen and McGovern [19] developed a model for theoretical 
decoupled Stirling cycle engine used for the analysis of the ideal 
Stirling cycle engine and its limits on its real-world realization. 

Thombare and Verma [20] provided a literature review of the 
important efforts carried out to develop the Stirling engine cycle 
and techniques used for engine analysis. They concluded that in or¬ 
der to achieve successful operation of the engine system with good 
efficiency, a careful design of heat exchangers with a proper selec¬ 
tion of the drive mechanism and the engine configuration are cru¬ 
cial. Their study indicated that a Stirling cycle engine working with 
a relatively low temperature with air as working fluid is potentially 
the attractive engines of the future, especially solar-powered low- 
temperature differential Stirling engines with vertical, double act¬ 
ing and gamma configuration. 

Karabulut et al. [21] conducted an experimental study of a Stir¬ 
ling engine with a lever-controlled displacer driving mechanism 
charged with helium. He concluded that the effect of charge pres¬ 
sure on the performance characteristics is positive up to a certain 
value and then negative. The relation between the hot end temper¬ 
ature and the power seems to be linear. For low and moderate heat 
sources including solar energy, the engine promises a reasonable 
performance. 

Kongtragool and Wongwises [22] designed and constructed a 
two single-acting, twin power piston and four power pistons with 
gamma-configuration and low temperature differential Stirling 
engine. They tested the engine with various heat inputs. The 
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variation of the engine torque, the shaft power and the brake ther¬ 
mal efficiency at various heat inputs with engine speed and perfor¬ 
mance were presented. They concluded that the engine 
performance increases while the heat input rises. The engine tor¬ 
que, the shaft power, the brake thermal efficiency, the speed, and 
the heater temperature also increase with the increasing heat in¬ 
put. Besides, the Beale number of this engine increases in parallel 
with the decreasing temperature ratio or with the increasing hea¬ 
ter temperature. 

Puech and Tishkova [23] conducted a theoretical investigation 
concerning the thermodynamic analysis of Stirling engine with lin¬ 
ear and sinusoidal variations of the volume. They demonstrated 
that the engine efficiency with perfect regeneration did not depend 
on the regenerator dead volume. In fact, this latter amplified the 
imperfect regeneration effect. 

Wongwises and Kongtragool [24] carried out a study on the 
power output determination of a gamma-configuration, low tem¬ 
perature differential Stirling engine. They demonstrated that the 
mean pressure power formula is most appropriate for the calcula¬ 
tion of a gamma-configuration, low temperature differential Stir¬ 
ling engine power output. Ross Yoke Stirling engine has many 
known advantages such as minimizing vibration and lateral forces 
acting on the pistons. Which lead to a more efficient and compact 
design with insignificant mechanical and thermal losses compared 
to other Stirling engines. In the present work, an attempt has been 
made to study this type of engine in order to improve its 
performance. 

According to the performed literature survey [25,15,26-30], 
Stirling engine of |3-type produces the highest indicated work 
whereas the y-type produces the lowest indicated work. Whereas, 
alpha Stirling engines drawn the minimum attention due to an old 
believe that it has low performance. So, there are few literatures 
concerning with the performance of this type of engines. However, 
it has many known advantages such as minimizing vibration and 
lateral forces acting on the pistons. Which lead to a more efficient 
and compact design with insignificant mechanical and thermal 
losses compared to other Stirling engines. In the present work, 
an attempt has been made to predict the performance behaviors 
of Ross-yoke Stirling engine in order to improve its performance. 

A detailed study of the influence of the Stirling engine design 
variables seems to be very useful. It will allow to optimize the dif¬ 
ferent engine components and to increase their performance. 
That’s why we developed a calculating code which is based on 
the dynamic model taking into account the losses in the different 
elements of the engine. This study demonstrated that a good opti¬ 
mization of the regenerator is necessary for this type of Stirling 
engine. 


2. Stirling engine classification 

2.1. Operating principle 

The engine integrated two volumes at variable temperatures re¬ 
lated together by a regenerative heat exchanger and auxiliary heat 
exchangers in order to realize the thermodynamic Stirling cycle. 

According to their geometric configuration, Stirling engines 
were basically classified into types. In fact, there are three config¬ 
uration systems: Alpha, Beta, and Gamma as shown in Fig. 1. 

These conditions described only the Stirling engine cylinder 
couplings which recognized the mode in which the displacer pis¬ 
ton and the power piston were linked, with respect to the connec¬ 
tion of the variable working spaces volume. These represent the 
spaces inside the engine cylinder where the working fluid is heated 
and cooled respectively in a closed volume. 


2.2. The different configurations of Stirling engines 

The majority of the literature concerning the Stirling engine was 
published sporadically since 1816 and there were no standardiza¬ 
tion or clarification of this terminology. This led to rely on diesel 
engine more than Stirling engine. A rich variety of literature en¬ 
gines is available, but no standard terminology concerning the best 
configuration exists. 

The Beta engine [31] includes the set of Stirling engine with a 
single cylinder pact and where the displacer and the power pis¬ 
ton are linked in tandem. The power is generated by the action 
of the pistons jointly. This type of configuration has sealing 
problems. Besides, it loses its advantages compared to multi-cyl¬ 
inder engines. 

The Gamma engine [32] is very comparable to beta-group since 
the power output is generated in the same way as in Beta engines. 
The only disparity is the fact that the two pistons move in separate 
cylinders. The considerable dead space that this configuration has 
represents the disadvantage of this group. 

The alpha engine [29] consists of two separate cylinders that 
each one has its own piston, either the displacer or the power pis¬ 
ton and three heat exchangers. Each heat source is linked to its 
own cylinder. The power output is produced by the separate mo¬ 
tion of the individual piston. 

The advantages of Alpha Stirling engine using the Ross yoke 
linkage (Fig. 2) are well known [33] because of the high power- 
to-volume ratio. However, the thermodynamic analysis of this en¬ 
gine was so restricted according to a literature overview since it is 
reported to necessitate high temperature but this problem was 
solved after the latest technological advances. 

Drive mechanisms cause a complicated problem for Stirling en¬ 
gines as discontinuous motion is necessary to achieve the volumet¬ 
ric changes that lead to a net power output [34]. 


3. Losses in a Stirling engine 

The energy loss in a Stirling engine is due to the thermodynamic 
processes and to the mechanical losses. The Compression and the 
expansion are not adiabatic. Load losses exist in the heat exchang¬ 
ers which are not ideal. 


3.1. Power loss by pressure drop SQ Pc h 


Pressure drops due to friction and to area changes in the heat 
exchangers can easily be modeled using steady flow correlations. 
These pressure drops lead to a dissipation of work within the heat 
exchangers leading to a reduction of the available work in the 
machine. 

The frictional drag force is expressed by [35]: 



where G is working gas mass flux (kg m” 2 s” 1 ), d is the hydraulic 
diameter, p is gas density (kg nrr 3 ), V is volume (m 3 ) and f f is the 
conventional fanning friction factor. 

A more suitable definition of the friction factor was 
proposed by Ureili [36], being so-called ‘Reynolds friction factor’, 
as follows: 

fr=f f Re (2) 


Re is the Reynolds number 


Re 


Gd 



The frictional drag force must be equal and opposite to the pres¬ 
sure drop force: 
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Fig. 1. Stirling engine configuration. 



rl, SQ Pchr 2 is the heat generated by pressure drop in the regenerator 
r2 and SQ Pchh is the heat generated by pressure drop in the heater. 

3.2. Energy Loss by internal conduction SQ Pc d 

Energy losses due to the internal thermal conduction between 
the hot parts and the cold parts of the engine through the solid ma¬ 
trix of the exchangers are taken into account. These losses are di¬ 
rectly proportional to the temperature difference between the 
two ends of the exchanger; they are given for the different 
exchangers [38]: 


dQpcdr = kcdry-(T r _ h - Tf_ r ) 

Cr 


dQpcdf — kcdf J (Jf-r T c ~f ) 

L f 


dQpcdh = Kdh T ~ (Jh-d — T r -h) 

Eh 


( 8 ) 


( 9 ) 


( 10 ) 


With SQ Pc d r , SQ Pcd f and SQpcdh representing respectively the conduc¬ 
tion loss in the regenerator, the cooler and the heater. 

The total conduction loss is then: 


F + ApA = 0 (4) 

Substituting Eqs. (2) and (3) into Eq. (1) we finally obtain: 

A 2 fr^ GV 

A P = ~ S 5 

Adp 

The internal heat generation which occurs when the gas is 
forced to flow against the frictional drag force, is given by [37]: 

= ( 6 ) 

m is the mass flow rate (kg s' 1 ). 

The total heat generated by pressure drop in the different 
exchangers is: 

dQpchT = dQ.Pchf "T dQpchrl ~T SQpchr2 T dQpchh (7) 

where SQ Pc hf is the heat generated by pressure drop in the cooler, 
SQpchr i is the heat generated by pressure drop in the regenerator 


dQpcdT — bQ-Pcdr "T bQpcdf "T dQ.Pcdh (11) 

To reduce the amount of energy lost from the system by inter¬ 
nal conduction, the kcdr, Ar, Lr and AT parameters should be 
optimized. 

3.3. Loss by external conduction 

The loss by external conduction is considered in the regenerator 
which is not adiabatic compared to the external side whereas for 
the other exchangers, the external losses are negligible. The energy 
stored by the regenerator at the time of the passage of gas from the 
space of expansion to the space of compression is not completely 
restored to this gas at the time of its return. So energy losses to¬ 
wards outside exist; they are determined from the regenerator 
effectiveness given by the following equation [39]: 

NUT 

~ 1 + NUT 

The energy lost by external conduction in the regenerator is so: 
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8Qpex = (}-E)SQr ( 12 ) 

With SQr = SQrt + <5Q.r2 

Qr! exchanged energy in the regenerator. 

The reduction of the energy lost by external conduction is very 
important since it represents the most significant loss in the Stir¬ 
ling engine. It depends on the regenerator effectiveness. 

3.4. Loss by shuttle effect SQpshti 

The losses by shuttle effect are due to the piston reciprocal ac¬ 
tion in the cylinder. The piston absorbs a heat quantity from the 
hot source and restores it to the cold source. Martini [37] proved 
that: 


dQpshtl 


0.4 Z 2 k pis D d 

JL d 


(T d ~ T c ) 



4. Design specification and concept 

4.1. Engine specification 

In order to avoid losses and to obtain high thermal efficiency 
the engine parameters should be optimized. The working fluid is 
air and the temperature difference between the heater and the 
cooler is about 600 °C only. 

The engine presented in Fig. 2 uses a Ross Yoke mechanism 
driving two pistons by the mean of yoke linkage [37]. The major 
feature of this is that the Ross Yoke is an ingenious mechanism 
for transferring dual piston motion into rotational motion. It has 
the advantage of minimizing lateral forces acting on the pistons 
leading to a more efficient and compact design. 


Table 1 

Volumes variations. 

Geometrical parameters 


Displacements 
Volume variations 


dV c 

d9 

dV e 

dO 


-A p r 

-A d r 


cos 9 + sin# 
cos# - sin# 



r sing cos 0 

be 

rsinOcosO 

be 


b i sin 0 = r cos # 

b9 = \Jb\(r — cos #) 2 

X = r sin # + b e 

Y c = r[sin # - cos #(b 2 /bi)] + b 6 
Y e = r[sin # + cos 0(jb 2 /fc»i)] + b e 

V c = V mc + A p (Y max -Y c ) 

V e = V me + A d (Y max -Y e ) 


Table 2 

Concepts and target performance. 


Parameters 

Values/type 

Engine type 

Alpha 

Working fluid 

Air 

Crank length 

r = 7.6 mm 

Yoke crank length 

hi = 29 mm 

Piston length 

b2 = 29 mm 

Displacement extremities 

Y m in = 17.75 mm 

Y max = 39.28 mm 

Swept volumes 

V bc - v mc = 24.42 cm 3 

Heat exchangers volumes 

V c = 5 cm 3 

V r = 6.75 cm 3 

14 = 4.78 cm 3 

Mean phase angle advance 

a = 90° 

Mass of gas in engine 

M= lg 

Hot space temperature 

T h = 923 K 

Cold space temperature 

T k = 350 K 

Frequency 

Freq = 41.72 Hz 


5. Dynamic model of Stirling engine 


4.2. Design concept 

The yoke drive mechanism does not generate sinusoidal volume 
variations. It is reported that the exact piston displacement func¬ 
tions are very intricate. 

The volume variations are derived from geometric consider¬ 
ations in Fig. 3 and Table 1. The major design concepts are listed 
in Table 2. 



There are many different ways to degrade the power produced 
by an ideal machine and in order to accurately predict the power 
and the efficiency; a clear understanding of the design compart¬ 
ments is required. 

A mathematical model takes into consideration the pressure 
drop in the heat exchangers [40-47]. 

Heat transfer and flow friction in the heat exchangers, i.e. the 
heater, the cooler and the regenerator, are evaluated using empir¬ 
ical equations under steady flow condition. 

No leakage is allowed either through the appendix gap or 
through the seals of the connecting rods. 

The gas temperature in the various engine compartments is 
variable. 

The cooler and the heater walls are maintained isothermally at 
temperatures Tpaf and Tpah. 

5.1. Dynamic model Equations 


The gas temperatures in the various 
lated from the perfect gas law: 

compartments are calcu- 

T Pc Vc 
c R ■ m c 

(14) 

T _ p f v f 
f Rm f 

(15) 

T P„V h 
" Rm h 

(16) 

T P d -V d 
d Rm d 

(17) 


Fig. 3. Geometric derivation of the Ross Yoke drive equation. 































































































154 


I. Tlili, S.A. Musmar/Energy Conversion and Management 68 (2013) 149-160 


The regenerator is divided into two cells rl and r2, each cell 
being associated with its respective mixed mean gas temperature 
T r i and T r2 expressed as follow: 


The work given by the cycle is: 


5W 

dt 



+ P d 


dVd 

dt 


T n = 

T r2 = 


P rl • V r i 

R ■ m r \ 

(18) 


Pr2 ■ Vr 2 

Rm r2 

(19) 

n 


The thermal efficiency given by the cycle is: 
W 

= Q^ 

The total engine volume is: 




An extrapolated linear curve is drawn through temperature val¬ 
ues T rl and T r2 defining the regenerator interface temperature T r _ ff 


T r _ r and T r _ h , as follows [32]: 


T 3T r i — T x2 

Tr - f = 2 

(20) 

T T r i + T r 2 

ir-r- 2 

(21) 

37,-2 — T r i 

T r - h = 2 

(22) 


According to the flow direction of the fluid, the interface’s tem¬ 
peratures: T c _ f , T f _ r , T r _ h and T h _ d are defined as follows [32]: 

If m c _ f > 0 then T c _ f = T c otherwise T c _ f = T f 


v T = v c + v f + v n + V r2 + v h + v d 

The other variables of the dynamic model are given from the en¬ 
ergy and mass conservation equation, applied to a generalized cell 
as follows: 

Energy conservation equation [33,34]: 

5Q + C p r £ m E -C p r s m s = P^ + C P ^^ (32) 

Since there is a varying pressure distribution throughout the 
engine, we have arbitrarily chosen the compression space pressure 
Pc as the baseline pressure. Thus at each increment of the solution, 
Pc will be evaluated from the relevant differential equation and the 
pressure distribution determined with respect to Pc. Thus it can be 
obtained from the following expression: 


If m f _ r > 0 then T f _ r = T f otherwise 7)_ r = T r _; 


If m r -h > 0 then T r _ h = T r _ h otherwise T r _ h = T h 

If m h _ d > 0 then T h _ d = T h otherwise T h _ d = T d 

where T c _ f is the temperature of the interface between compression 
space and the cooler, T f _ r is the temperature of the interface be¬ 
tween the cooler and regenerator, T r _ h is the temperature of the 
interface between the regenerator and the heater, T h _ d is the tem¬ 
perature of the interface between the heater and the expansion 
space. 

The matrix temperatures are thus given by: 


„ „ A Pf 

P { = P C + / 

(33) 

P rl =P / + (A ^ + 2 APrl) 

(34) 

P r 2-P r ,+ (APrl 2 APr2) 

(35) 

P h = P, + (APrt j ^ 

(36) 

P d - P„ + APh 

(37) 


dTpa n _ SQ n 
dt ~ c pr dt 



Applying energy conservation equation to the various engine 
cells, we obtain: 


dTpa r2 _ SQ r2 
dt c pr dt 

where Cpr is the heat capacity of each cell matrix (J K -1 ), Qrl is the 
quantity of heat exchanged to the regenerator rl (J), Qr2 is the 
quantity of heat exchanged to the regenerator r2 (J), T Pan is the ma¬ 
trix temperatures in the regenerator rl (K) and T Par2 is the matrix 
temperatures in the regenerator r2 (K). 

By taking into account the conduction loss in the exchangers 
and the regenerator effectiveness, the power exchanged in the var¬ 
ious exchangers is written: 


8Qf = hfApaf(Tpaf — Tf) 

(25) 

^Q.r2 — d r 2^par2 (Tpar2 ~ Tr2) 

(26) 

^Q.rl — d r ^Apar1 (Tpari ~ Trl) 

(27) 

SQh = h h Ap ah (T pah -T h ) 

(28) 

The heat transfer coefficient of exchanges h f , h n , 
only available empirically [29]. 

The total exchanged heat is: 

h r2 and h h are 

<5Q. = SQf + ^Q.rl T" 3Qr2 ~h &Q.h 

(29) 


CpTc-fifics — r ^Q>P ( 


^ + c,v, dp - 


dt 


dt 


(5Q/ - SQpchf + C p T c _firif E - C p Tf_ r rhfs = 


C v V f dP c 


R dt 


^Q.r 1 ^Qpchr 1 "T CpTj_ r rh r i£ CpT r _ r 17lrlS — 


CyVr 1 dP c 
R dt 


^Q.r2 dQ.Pchr2 H - CpT r _ r /?l r 2£ CpT r -h1Tlr2S 


C v V r 2 dP c 


R dt 


&Qh ^Q.Pchh H - CpTr-hdlhE Cp^h-e^hS 


CyV h dP c 


R dt 


rrr . 1 (r n&V d r Xf dP c \ 

C P T h . d m d - R (C p P d d( , + C v V d d[ j 

Summing Eqs. (39)-(44) we obtain the pressure variation: 


dPc 

dt 


1 

cvw 


R(SQ — SQpchr ) - 


SW 

dt 


Mass conservation equation: 


(38) 

(39) 

(40) 

(41) 

(42) 

(43) 

(44) 


M = m d + m c + m/ + m r + m h 
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Table 3 

Calculations of error analysis. 


Cycle iteration 

1 

2 

3 

4 

5 

6 

7 

8 

Power by cycle iteration 

23.8807 

22.5903 

22.2088 

22.1731 

22.1701 

22.1698 

22.1698 

22.1698 

Efficiency by cycle iteration 

0.3760 

0.3030 

0.2949 

0.2942 

0.2942 

0.2942 

0.2942 

0.2942 

Error for power, e p 

1.2904 

0.3815 

0.0387 

0.003 

0.003 

0.0000 

0.0000 

0.0000 

Error for efficiency, s e 

0.0730 

0.3428 

0.0081 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


The mass flow in the various engine compartments are given 
from the conservation equations of energy (39)-(44): 
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5.2. Method of solution 
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The Stirling engine tested has a Ross Yoke linkage system, as 
shown in Fig. 2. The geometrical parameters of this engine are gi¬ 
ven in Table 1. 

The operating conditions are as follows: working gas: air at a 
mean pressure of 2.195 bar; frequency: 41.72 Hz; hot space tem¬ 
perature: T Pa p = 923 I<; cold space temperature: T Paf = 350 K. 

The dynamic model developed in this work including thermal 
irreversibilities was applied and tested on the Ross Yoke Stirling 
engine. 

The independent differential equations obtained in paragraph 4, 
are solved simultaneously for the variables P c , m c , T r i, W, etc. 

The problem thus reduces to the simultaneous solution of a set 
of ordinary differential equations. The simplest approach to solve 
this set of ordinary differential equations is to formulate it as an 
initial-value problem, where the initial values of all the variables 
are known and the equations are integrated from this initial state 
over a complete cycle, bringing the pistons back to their initial 
positions. 

The corresponding set of differential equations is expressed as 
dYldt = F(t,Y). The objective is to find the unknown function Y(t) 
which satisfies both the differential equations and the initial con¬ 
ditions. By assigning arbitrary initial conditions for the seven vari¬ 
ables to be integrated and integrating the equations through 
several complete cycles, the cyclic steady state can be attained 
when the respective values at the beginning of the cycle and at 
the end of the cycle are equal. The residual regenerator heat Qr 
at the end of the cycle, should be zero. The system of equations 
is solved numerically using the classical fourth order Runge-Kutta 
method, cycle after cycle until periodic conditions are reached. The 
calculations of error analysis are cited in Table 3. 

To validate the numerical method used in the computation, the 
results are compared with those obtained by Urieli and Berchowitz 
[36] for the same conditions (adiabatic models) of the GPU-3 en¬ 
gine data. The comparison shows a good agreement, as shown in 
Fig. 4. 


- Btrchowitx rt tuh« (Adiabatic modal) 

XXX Our result5 (Adah*!* model) 

Fig. 4. Validation of the computational method. 



6. Results and discussion 

Stirling engines are mechanical devices theoretically based on 
the Stirling cycle. That is a highly idealized thermodynamic cycle. 
It consists of four processes which combine to form a closed cycle: 
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two isothermal and two isochoric processes. The surface of the PV 
diagram is equal to the network produced by the engine. 

The maximal pressure reached in the engine is of 2.507 bars at 
the end of the compression phase, and the mean pressure is 2.195 
bars. 

The surface of the PV diagram shown in Fig. 5 is equal to the 
network produced by the engine. 

We notice that the work calculated by the dynamic model is 
lower than that found by the isothermal model since the PV dia¬ 
gram surface is narrowed. This figure shows also the effect of the 
sinusoidal movement of the displacer piston on P-V diagram. It 
is noticed that the average pressure is increased compared to the 
adiabatic case, due to the effect of temperature variation. 

The regenerator is one of the most complicated components of 
the Stirling engine and is the major limiting factor in Stirling en¬ 
gine performance. It acts as a thermal sponge alternatively accept¬ 
ing heat from working fluid as it passing through it toward the 
cooler segment and rejecting heat to the working fluid in its way 
back. 

6.2. Heat energy and mass flow in different heat exchangers 

Heat exchangers are key components in the Stirling engines 
which are the heater, regenerator and cooler which exchange heat 
to and from the engine. The study of the heat transfer behavior 
within and across the Stirling engine may provide an adequate pre¬ 
diction for overall system thermal efficiency in practice. 

Figs. 6 and 7 show the heat energy exchanged and torque vari¬ 
ation, according to the dynamic model. Regenerators are typically 
made up of porous structures, usually referred to as the packing 
material, which may lead to complex flow and heat transfer behav¬ 
iors of the working fluid while passing through the regenerator. 
The heat transferred in the regenerator is very significant almost 
ten times than that of the net work done per cycle. The heat trans¬ 
fer coefficient is much lower for gases than for liquids, thus the 
vast surface area in a regenerator greatly increases heat transfer. 
It is concluded that the engine performance depends critically on 
the regenerator effectiveness and its ability to accommodate the 
high heat flow. 



0 50 100 150 200 250 300 350 

Cycle angle (degrees) 

Fig. 6. Heat energy variation. 



Cycle angle (degrees) 

Fig. 7. Torque variation. 

When the cycle angle increase the amount of heat exchanged 
increases until it reaches a maximum equal to 24 J. A regenerator 
is used to increase the efficiency of a Stirling engine. 

Fig. 6 prove also that the regenerator is well designed since the 
energy rejected by the gas to the regenerator matrix in the first half 
of the cycle is equal to the energy absorbed by the gas from the ma¬ 
trix in the second half of the cycle. Therefore, High theoretical effi¬ 
ciency associated with Stirling engine compared with internal 
combustion engine is mainly due to the absence of intake and ex¬ 
haust process (i.e. no energy exhaust from the Stirling engine). 

The power and the thermal efficiency calculated by this model 
are closer to the real power and thermal efficiency of the prototype 
than other model. Since the average power is 48.79 W and the en¬ 
gine effectiveness resulted is up to 41%. 

It is noted that the flow in the cooler reaches a value greater 
than the flow in the regenerator which reaches a higher value than 
the heater. This is due to the location of the cooler that is close to 
the engine piston. 

6.2. Pressure drop and energy dissipation in the different heat 
exchangers 

Stirling engines have pressure losses at heater, regenerator, 
cooler and their connection parts. The design and configuration 
of these heat exchangers affect the engine performance. We Notice 
from Fig. 8 showing the variation of pressure drop in the heater, 
the regenerator and the cooler that the regenerator is responsible 
of the undesirable pressure decrease, which attains 0.064 bars be¬ 
cause of its geometrical and physical aspects. Whereas the pres¬ 
sure drop in the other exchangers reaches a maximum 
0.00268 bar. Large effective heat transfer area and few dead vol¬ 
umes cause high pressure drop across the regenerator. This in¬ 
crease of pressure drop across the regenerator reduces the 
output power. 

Also from Fig. 9 we find that the energy dissipation in the regen¬ 
erator is so high compared to the other exchangers. 

It is also noticed that approximately 51.24% of power is lost by 
the regenerator compared to the net power of the engine. However 
its rule is major for the engine output. 

6.3. Optimization of the Stirling engine performances 

Stirling engines are known as high efficiency engines. In prac¬ 
tice due to some reasons their efficiency reduces. Some techniques 
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Cycle angle (degrees) 

Fig. 8. Pressure drop variation. 


Figs. 10 and 11 show the effect of the heater temperature on the 
engine power and efficiency for different frequencies. It is noted 
that the engine performances increase according to the heater 
temperature and the power increases with the operating 
frequency. 

We note that when Tpah < Tpaf\ the machine can function as a 
refrigerator; when the expansion temperature decreases, it is then 
necessary to increase the mechanical power. We demonstrated 
that we can assemble the machines in two ways: one like a refrig¬ 
erator and the other like an engine. 

The curves indicate that the increase of power with heat source 
temperature is limitless. Therefore, the progress made in the mate¬ 
rial technology especially the resistance to high temperature, will 
improve the performance of Stirling engines. 

We notice also that the increase of the heater temperature in¬ 
duces a continuously increase in the heat exchanged. 

Although the engine losses increase when the heater tempera¬ 
ture rises [13], the performances of the engine augment, this is ex¬ 
plained by the increase of the exchanged energy between the 
matrix and the working fluid of the regenerator. The temperature 
of working fluid entering the regenerator is not constant because 
the pressure, density and velocity of working fluid vary over a wide 



Cycle angles (degrees) 


Fig. 9. Dissipation flow in different heat exchangers. 



Heater temperature (K) 

Fig. 10. Temperature effect on engine power. 


can give an insight into the detailed interaction of factors deter¬ 
mining performance, and can be used to produce a significant 
improvement in the predicted performances of this particular en¬ 
gine. One of the most important sources of efficiency reductions 
is imperfect heat transfer. The search for an engine cycle with high 
efficiency has led to improve performance analysis and optimiza¬ 
tion of Stirling engine which depends on the losses reduction and 
the determination of the adequate geometrical and physical 
parameters. In fact, Stirling engine loses efficiency due to large dif¬ 
ferences in the temperature of the working fluid and the heating 
and cooling spaces. Usually the design point of a Stirling engine 
will be somewhere between the two limits of: maximum efficiency 
point; and maximum power point. Therefore, it is required that for 
successful operation of engine system with good efficiency a care¬ 
ful design of heat exchangers, proper selection of drive mechanism 
and engine configuration is essential. 


6.3.2. Effect of the heater temperature 

The heater transfers heat from external source to the working 
fluid contained within the engine working space. 



Heater temperature (K) 


Fig. 11. Temperature effect on engine efficiency. 
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range. The effectiveness of the regeneration process largely de¬ 
pends upon the thermal capacity of the material also. 

6.3.2. Effect of the fluid mass 

Figs. 12 and 13 show the fluid mass effect on the engine power 
and efficiency. It is noted that the power increases according to the 
mass and that the efficiency reaches a maximum. When the mass 
increases, the reduction in the output is due to the load loss 
increase. 

An increase of the total mass of gas in the engine leads to a rise 
in the density, the mass flow, the gas velocity, the load, and pres¬ 
sure function. Therefore, this increase in the total mass of gas in 
the engine will lead to a more energy loss by pressure drop [16]; 
however, the engine power increases and the efficiency reaches a 
maximum value. 

The engine efficiency decrease occurring when the mass of fluid 
increases beyond a certain value is a consequence of the reduction 
of the maximum temperature of the working fluid secondary to the 
transfer capacity limitation in the heater and the regenerator. 
When the mass increases further, the decrease of the efficiency is 
due to an increase of the pressure loss and the limitation of the 
heat exchange capacity in the regenerator and the heater. 


The two Figs. 12 and 13 show that more we increase the tem¬ 
perature more the engine performance increases and that the opti¬ 
mal value of the fluid mass decreases when temperature increases. 
It is known that the whole working fluid mass is divided between 
the cold and hot spaces. The mass also takes part in every process 
consistently and the process makes up a cycle that includes com¬ 
pression and expansion. 

6.3.3. Effect of the swept volume ratio of compression and expansion 
on the power engine 

Figs. 14 and 15 show the effect of the swept volume ratio of 
compression and expansion (t = Vc/Vd) on the power and the effi¬ 
ciency. The curves show clearly that there is a well defined optimal 
value of t at which the efficiency is maximum. The comparison of 
the two curves for ( Tpah/Tpaf) = 0.25 and ( Tpah/Tpaf') = 0.33 shows 
that the optimal value is not constant but changes from 0.8 when 
(' Tpah/Tpaf) = 0.25-1.1 when ( Tpah/Tpaf) = 0.33. The changes of the 
dead volume ratio produce also an adjustment of the optimal value 
of t. Thus, there is not a better fixed value t. The decrease of the 
output when t increases beyond a certain value is probably 
compared to the instantaneous distribution of the fluid masses in 
the different exchangers during the cycle. A better mechanical 
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Fig. 12. Effect of working fluid mass on engine power. 


Fig. 14. Effect of the swept volume ratio on the engine power. 
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Fig. 13. Effect of working fluid mass on engine efficiency. 


Fig. 15. Effect of the swept volume ratio on the engine efficiency. 
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efficiency can be obtained by taking a swept volume ratio smaller 
than an optimal value which vary according to the ( Tpah/Tpaf) 
while resulting in only a relatively small reduction in shaft output. 


Increase in dead volume within the Stirling engine system re¬ 
sults in a loss of power but not necessarily a reduction in efficiency. 


6.3.4. Effect of the dead volume ratio 

Figs. 16 and 17 show the effect of the regenerator volume rate 
compared to the heater (cp = V r IV h ) on the performances of the en¬ 
gine. In a Stirling engine, the volume of the heater and regenerator 
are dependent and function of its power. These performances are 
function of the temperature ratio Tpaf/Tpaf. 

It is very clear that the increase of dead space above the mini¬ 
mum required reduces the power parameter. Dead space in the 
heater would be reduced but the temperature drop between source 
and working fluid would be increased in inverse ratio. Further in¬ 
crease in efficiency can be obtained by maintaining the dead space 
as small as possible. Thus, the regenerator volume should not be 
small in order to ensure a good thermal change and it should not 
be large so that it does not represent a dead volume which penal¬ 
izes the output of engine. Therefore, engine power increasing stea¬ 
dily as regenerator length is reduced. This is to be expected; power 
would be maximized by doing away with the regenerator alto¬ 
gether, provided enough heat could be supplied to make up for 
the resulting inefficiency. 
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Fig. 16. Effect of the dead volume ratio on the engine power. 
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Fig. 17. Effect of the dead volume ratio on the engine efficiency. 


7. Conclusion 

According to the performed literature survey, alpha Stirling en¬ 
gines using the Ross Yoke linkage drawn the minimum attention 
due to an old judgment that it has low experimental performance. 
However, it has many known advantages such as minimizing 
vibration and lateral forces acting on the pistons. Which lead to a 
more efficient and compact design with insignificant mechanical 
and thermal losses compared to other Stirling engines. In the pres¬ 
ent work, an attempt has been made to study this type of engine in 
order to improve its performance. Since this paper represents an 
application of the thermodynamic modeling in order to optimize 
the operating conditions of the Stirling engine, so that we obtain 
the best mode of power and performance. The parametric study 
of the engine performances highlighted the importance of the 
choice of some parameters (Tpah, M, t, cp, co) and their influence 
on the engine performance indices (Power, efficiency). We ana¬ 
lyzed also the regenerator volume effect on the engine perfor¬ 
mance. The results obtained on the basis of the total model 
taking into account the thermal and mechanical losses show the 
fact that the regenerator parameters (thermal effectiveness, vol¬ 
ume, losses) have a great impact on the performances of the Stir¬ 
ling engine. All these results are in agreement with the known 
experimental data for the Stirling engine [18,24,32,29,45,47]. 

The engine performance can be improved by increasing the pre¬ 
cision of the engine parts and the heat source efficiency. The engine 
performance depends critically on the regenerator effectiveness 
and its ability to accommodate the high heat flow; it should be in¬ 
creased if a better working fluid e.g. helium or hydrogen is used in¬ 
stead of air. 
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